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Abstract. The fraction of AGN with photoelectric absorption in the X-rays ranging from of 10^^ up to about lO"'* cm“^ 
(Compton-thin) appears observationally to be anticorrelated to their luminosity G-. This recently found evidence is used to 
investigate the location of the absorbing gas. The molecular torus invoked in the unified picture of AGN, while it can be 
regarded as confirmed on several grounds to explain the Compton-thick objects, do not conform to this new constraint, at least 
in its physical models as developed so far. In the frame of observationally based evidence that in Compton-thin sources the 
absorbing gas might be located far away from the X-ray source, it is shown that the gravitational effects of the black hole 
(BH) on the molecular gas in a disk, within 25-450 pc (depending on the BH mass, from 10® to 10® Mq), leads naturally to 
the observed anticorrelation, under the assumption of a statistical correlation between the BH mass and L^. Its normalization is 
also reproduced provided that the surface density, E, of this gas is larger than about 150-200 Mq pc“^, and assuming that the 
bolometric luminosity is one tenth of the Eddington limit. Interestingly, the required values are consistent with the value of the 
300 pc molecular disk in our own galaxy, namely 500 Mq pc“^. In a sample of nearby galaxies from the BIMA SONG survey, it 
is found that half of the objects have central S larger than 150 Mq pc“^. Given the simplicity of the proposed model, this finding 
is very encouraging, waiting for future higher resolution surveys in CO on more distant galaxies. 
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1. Introduction 

Similar to dust extinction in the optical, photoelectric absorp¬ 
tion in X-rays provides clues on the environment of Active 
Galactic Nuclei (AGN). In general the latter provides informa¬ 
tion which is more straightforward to interpret than the former, 
because it depends a) on atomic properties and consequently 
only on the chemical composition of the medium, b) on the 
adoption of a shape for the X-ray continuum, which is empir¬ 
ically known to follow a simple Power Law (PL). The values 
of the spectral photon index, F, display a rather modest disper¬ 
sion (e.g. Perola et al. 2002; Piconcelli et al. 2005); further¬ 
more, in any single object under consideration, and provided 
that the absorption is Compton-thin (i.e. with Nh not exceed¬ 
ing cr"*=L5xl0“^^ cm“^), the exact value of F can be directly 
measured, given an adequately wide observational band and 
sufficient statistics for the spectral counts. Here, though, it is 
more important to stress the limitations, which are mainly of 
three sorts: 

1. The source of the PL photons is, for most practical pur¬ 
poses, point-like, and therefore it remains questionable how 
to relate the total absorbing matter density (the absorbing 
column Nff) to bi-dimensional information (either direct or 
indirect) on the circumnuclear matter from observations in 
other wave-bands. In other words, usually it cannot be ex¬ 
cluded that the observed column might be pertaining to that 


unique line of sight, rather than to the overall configuration 
of the circumnuclear matter. 

2. While it is spectroscopically possible (as amply demon¬ 
strated even before the use of very high resolution instru¬ 
ments) to disentangle the imprints of a multi-phase (in ion¬ 
ization terms) absorber, it is by no means immediately pos¬ 
sible to recover the space distribution of the phases along 
the line of sight. When the ionization is attributed to the 
photoelectric effects of the PL source as a reasonable hy¬ 
pothesis (better still when it can be observationally proven 
not to be of collisional origin), it is generally regarded as 
most likely that the high ionization gas is very close to the 
centre. This leaves however quite open the possibility that 
the low (in practice, indistinguishable from neutral) ioniza¬ 
tion gas might lie either much further away or equally close 
but dense enough to remain neutral. 

3. When N// is larger than about 10^"^ cm“^, the column 
is Compton-thick, that is a diffusive process adds to the 
absorption process to deplete the specific intensity along 
the line of sight. A truly thick absorber (t^/, » 1) re¬ 
duces the intensity by several orders of magnitude all the 
way throughout the hard X-ray band, and renders in fact 
its quantitative evaluation impossible, except for a “lower 
limit” (e.g. Matt et al. 1999). 

It must be emphasized that, in general, AGN spectra con¬ 
tain another component in addition to the PL, which becomes 

evident beyond 8-10 keV. It is due to reflection from Compton- 
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thick, neutral gas, and, as expected on physical grounds, it is 
accompanied by a strong fluorescent iron line. From a pure ob¬ 
servational point of view, this component is often the only vis¬ 
ible one below ~ 10 keV when the PL is photoelectrically ab¬ 
sorbed by N// of about 10^“* cm“^ (corresponding to a turnover 
in energy at about 10 keV), or greater. Its intensity is typically 
a few percent of the primary one in the 2-10 keV band. In such 
systems, it is regarded as very reasonable to assume that the 
absorbing and the reflecting gas belong to one and the same 
circumnuclear structure. Historically, it was the spectropolari- 
metric measurements in the optical band of NGC 1068 which 
led Antonucci & Miller (1985) to propose the so-called “uni¬ 
fied model”. This model (Antonucci 1993) envisages the exis¬ 
tence of a thick and dusty “torus” of molecular gas around the 
nucleus, extending out to a few tens of pc at most (Risaliti et 
al. 1999) and with a substantial covering factor to account for 
the large ratio of type 2 to type 1 Seyferts, at least in the local 
Universe (4:1, if also Seyfert 1.8 and 1.9 are included; Maiolino 
& Rieke 1995; Ho et al. 1997; it is worth recalling that the lo¬ 
cal Universe is dominated by relatively low luminosity objects 
as compared to the Quasi Stellar Objects, QSO). In its basic 
form, the model has been widely confirmed, most notably by 
measurements in the hard X-ray band, which showed that most 
type 2 Seyferts are X-ray absorbed. The situation, however, is 
likely to be more complex than the zero-th order model would 
predict. On one hand, there is increasing evidence of a sub¬ 
stantial population of optically inactive, obscured AGN, both 
in the local (Maiolino et al. 2003) and more distant (Comastri 
et al. 2002; Brandt & Hasinger 2005) Universe. This implies 
that the ratio bewteen X-ray obscured to unobscured AGN may 
be larger than that between type 2 and type 1 Seyferts. On the 
other hand, other absorbing regions, besides the torus, may ex¬ 
ist in the complex environment of AGN: dust lanes (Malkan 
et al. 1998), starburst regions (Weaver 2002), the galactic disk 
(Maiolino & Rieke 1995). Indeed, it has been suggested (see 
e.g. Matt 2000, 2004 and references therein) that Compton- 
thick absorption (which, in the local Universe, is observed in 
about half of optically selected type 2 Seyferts, Risaliti et al. 
1999, Guainazzi et al. 2005, the percentage rising when “elu¬ 
sive” AGN are also taken into account, Maiolino et al. 2003) 
is related to the torus, while Compton-thin absorption is due 
to one (or more) of the other possible absorbing regions. (If 
a dust/gas ratio typical of the cold interstellar medium of our 
Galaxy is adopted, Compton-thin absorption, i.e. N// in the 
1022-10^4 cm“2 range, is sufficient to explain the extinction of 
the BLR lines). 

This paper will address the question of the location of the 
absorber in the light of new results from X-ray surveys in the 
2-10 keV band. The sensitivity reached in this band with the 
XMM-Newton and Chandra satellites made it, for the first time, 
possible to study AGN of both type 1 and 2, with sufficient 
statistics to investigate the N// distribution as a function of the 
intrinsic (PL) luminosity and of the cosmic epoch out to z about 
4. 

Surveys limited to 10 keV on the high energy side (and 
with instrument sensitivity dropping rapidly beyond 5-7 keV) 
are doomed to miss almost completely the Compton-thick ob¬ 
jects. Thus the results must be read as pertinent, among the 


AGN classified through optical spectroscopy as type 2, only 
to those with Nh less than 10^"^ cm 2. Basically, the question 
addressed in this paper is whether the new results can be un¬ 
derstood in the context of the unified model based on just one 
element, the torus mentioned above, or, if not, if one can con¬ 
ceive of just one realistic (astrophysically) additional element 
able to provide a positive answer, with the least possible num¬ 
ber of assumptions. In particular, the new contribution of this 
paper is the use of one further constraint, namely the decrease 
of the fraction of Compton-thin obscured AGN with increas¬ 
ing X-ray luminosity, an anticorrelation which have recently 
emerged from the statistical analysis of the abovementioned 
types of survey (Ueda et al. 2003; La Franca et al. 2005). 

In Sect. 2 this new constraint is described. In Sect. 3 the cur¬ 
rent physical models of the torus are discussed and shown to be 
unable to account, even qualitatively, for the anticorrelation. In 
Sect. 4 it is shown that, if the obscuration is due to molecular 
gas in galactic disks, then the gravitational effects of the black 
hole on its spatial distribution within 25^50 pc (depending on 
the mass) leads naturally to an anticorrelation. An encourag¬ 
ingly reasonable quantitative match is achieved under two con¬ 
ditions: that the molecular mass surface density exceeds ~150 
Mq pc“2 within 25-450 pc, and that in a statistical sense the 
bolometric luminosity is proportional to the black hole mass, 
and equal to about 10% of the Eddington luminosity. In Sect. 5 
these two conditions are discussed in the light of the evidence 
available, and a positive conclusion is drawn. 

2. A new observational constraint: the 
luminosity-absorption anticorrelation 

A statistical investigation on the photoelectric absorption prop¬ 
erties of Compton-thin AGN was first obtained, down to a flux 
limit F(2-10 keV) of about 10“'^ erg cm“2 (which corre¬ 
sponds to about 25% of the X Ray Background, XRB) with 
imaging instruments operating up to 10 keV onboard ASCA 
and BeppoSAX (e.g. Ueda et al. 1999; Fiore et al. 1999). The 
main problem with these pilot surveys is the error boxes at 
fluxes close to the limit, whose size left substantial uncertain¬ 
ties on the optical identification, especially for the type 2 AGN 
(see La Franca et al., 2002, for a discussion on this point). 
The situation improved greatly in the last four years thanks to 
Chandra and XMM-Newton, both of which image the sky up 
to 10 keV with much higher sensitivity and much better angu¬ 
lar resolution (at the arcsecond level in the case of Chandra). 
The combination of these two qualities have led to positionally 
based identifications, followed by spectroscopic identifications 
of the counterparts, which are regarded as fairly secure down 
to magnitudes R about 24. The X-ray flux limit achieved with 
a high fraction of identified sources is approximately 10 erg 
cm"2 (e.g. La Franca et al. 2005). With the instrument sen¬ 
sitivity covering the spectra down to 0.5 keV, along with the 
knowledge of z, it is possible to estimate both the intrinsic value 
of N// and the X-ray luminosity L^ corrected for absorption. 
After accumulating a sufficiently large sample to cover the L^-z 
plane, one can investigate simultaneously the X-ray luminosity 
function of AGN (irrespective of their being optically classified 
as type 1 or 2), its cosmic evolution, and the distribution of N// 
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as a function of and z. In the following, will be referred 
to as the intrinsic luminosity in the 2-10 keV band. As noted 
in Sect. 1, truly Compton-thick sources remain basically un¬ 
sampled, thus the constraint emerging from these studies apply 
only to sources with N jf ranging from (nominally) zero to about 
10^"^ cm“^. With the statistics achieved so far it appears prema¬ 
ture to observationally quantify the functional form of the N// 
distribution. At present it is convenient to adopt a conservative 
approach, by using the ratio ^ = (number of sources with log 
(Nh) > 22)/(total number of sources). Ueda et al. (2004) have 
shown that ^ is dependent on that is ^ = ^(L^), where ^ de¬ 
creases with increasing L^. Moreover, the same authors reached 
the conclusion that this property appears to be quantitatively in¬ 
dependent from z. Recently, with a larger sample, which over¬ 
laps the one utilized by Ueda et al. (2004), but yields a wider 
coverage of the plane. La Franca et al. (2005) confirmed 
the first finding, but not the second. Namely they find evidence 
of a z dependence, thus formally ^(L^, z), where ^ increases 
and becomes shallower as z increases. 

In this paper only the dependence of f on in the local 
universe, as obtained by La Franca et al. (2005) from z=0 to 
z= 0 . 1 , will be used as the novel constraint, while the depen¬ 
dence on z will be briefly commented upon and left as the sub¬ 
ject of further investigations. 

The next Section is devoted to investigate whether this an¬ 
ticorrelation can be understood in the frame of physical models 
of the torus. 

3. The case of the “thick torus” 

Several authors (Krolik & Begelman 1986, 1988, Beckert & 
Duschl 2004, Vollmer, Beckert & Duschl 2004, Beckert et al. 
2004) have addressed in physical terms the dynamical and ge¬ 
ometrical structure of the ’’thick torus”, which is assumed to 
be cylindrically symmetric. If R is the distance from the AGN 
centre (the site of the black hole, BH, with mass Mbh), and 
H(R) is the scale height (vertical geometrical thickness), these 
models are often depicted (see Fig. 1) as if the ’’torus” were 
sharply confined within 2H{R), with all the obscuring matter 
within and nothing outside. This type of descriptive picture ap¬ 
plies well to the separation of the Compton-thick sources from 
the rest. One could imagine a smoother transition of the obscur¬ 
ing matter density, in order to accomodate also the absorbed, 
but Compton-thin sources. The question is then the following: 
do these models, on the grounds of their physical principles 
and self-consistency, behave in a way that, without additional 
ad hoc assumptions, comply with the constraints described in 
Sect. 2? In these models H is an increasing function of R, sim¬ 
ply as a consequence of the gravity exerted by the BH, in addi¬ 
tion to that exerted by a concentric and, in first approximation, 
spherically symmetric stellar structure. In order to ensure the 
required scaleheight together with the most effective obscura¬ 
tion, the matter is assumed to be in the form of dense molecular 
clouds, endowed with the adequately large velocity dispersion. 
The clouds are considered to be either self-gravitating (Beckert 
& Duschl 2004, Vollmer, Beckert & Duschl 2004) or in pres¬ 
sure equilibrium with a hot phase. The clouds can closely inter¬ 
act with each other or be subject to tidal shear, both effects lead¬ 


Type 1 



Fig. 1. Geometry of the ’’thick torus” (Nenkova et al. 2002) 


ing to their destruction; they set stringent limits on the cloud 
size (hence their internal density) as well as on their number 
density as a function of R. Beckert et al. (2004) have shown 
that the geometrical thickness of the torus is given by: 


H = 


2R^M 

CsM{R) 


( 1 ) 


where M is the total mass accretion rate, M{R) the mass within 
the radius R and Cs is the sound speed. 

In order to maintain the obscuration effective as required by 
observations (see Sect. 1), M must be much larger than MEdd, 
namely most of the matter contributing to M must be lost, and 
only a fraction of it {Mium) is allowed to feed the AGN power 
(see Fig. 4 in Beckert & Duschl 2004.). If increases with 
M, which seems reasonable to assume, H/R, hence the covered 
fraction of Att, would increase with the luminosity, contrary to 
the novel constraint. 

One might of course read this conclusion as suggesting that 
something, to be physically demonstrated, is acting in the sys¬ 
tem to ensure anticorrelation between M and Mium- So far no 
account has been taken of the feedback on the state of the mat¬ 
ter due to the irradiation by the AGN itself. The ionization- 
evaporation effects set a limit on the minimum distance to 
which an absorbing cloud can survive. A useful quantification 
follows (Krolik & Begelman 1988) from equating the “evapo¬ 
ration time-scale” and the “accretion time-scale”: 


R. 


0.79Mgfls^5L44 

~ 


pc 


( 2 ) 


where L 44 is the central ionizing luminosity in units of 10“*^ erg 
s“', Mgas ,5 is the gas mass in units of IQ^Mq, Ncip,A is the col¬ 
umn density of a single cloud in units of 10 ^“* cmr'^, 7$ is the 
temperature at the sonic point of the evaporating flow in units 
of 10^ K, and Rout is the outer radius of the configuration. In a 
purely geometrical model of the “torus”, described as a cylin¬ 
der with a cavity (Lawrence 1991), the radius of the latter is 
bound to increase with the luminosity. Consequently, the cov¬ 
ered fraction of Att would decrease. However this would be an 
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acceptable answer to the anticorrelation described in Sect.2 if 
one were able to demonstrate that physically HIR can be main¬ 
tained constant (or decrease) versus R. The physical models 
summarized above, if M(R) increases with R less than linearly 
(as it is the case when this term is dominated by the BH mass) 
show instead that the pure effect of the BH gravity involves a 
concave geometry (see Fig. 3 in Beckert & Duschl 2004), hence 
the covered fraction is basically determined by H at Rout- This 
implies that the fraction of obscured sources would be inde¬ 
pendent of the luminosity up to a certain value at which the 
illumination is so strong to forbid the formation of the torus. 
This is thence an on-off effect, which could be relevant to the 
issue of the fraction of Compton-thick sources as a function of 
their luminosity: this issue at present lacks an observationally 
secure point of reference, and will not be further discussed in 
this paper. 

4. Obscuration by molecular gas in the galactic 
disk 

On the basis of a statistical analysis of the correlation between 
the spectroscopic optical classification of nearby AGN and the 
inclination to the line of sight of the host (typically spiral) 
galaxies, Maiolino & Rieke (1995) found the following. For 
the extreme type 2 there seems to be no correlation. On the 
contrary, types l-r 1.24-1.5 are mostly found in face-on galax¬ 
ies, and the intermediate types, 1.84-1.9, are mostly found in 
edge-on galaxies. They therefore suggested that the obscura¬ 
tion in the intermediate AGN types is associated with the gas 
in the galactic disks. Risaliti et al. (1999) showed that these 
types, when observed in X-rays, correspond to Compton-thin 
AGN. 

It remains to be seen if the anticorrelation discussed in Sect. 
2. could be explained within this scheme. The gas content at 
large in the disk of spiral galaxies correlates with their mor¬ 
phological type. In particular, the mean HI surface density over 
the disk out to the same isophotal radius is about twice in Sc 
that in Sa galaxies (Roberts & Haynes, 1994). However, the 
concentration of the gas is generally associated with the spiral 
arms. So, it should increase towards the center in the Sa type, 
where the arms are more tightly wound, and decrease in later 
types. Furthermore, since the far infrared emission (FIR) ap¬ 
pears to correlate with the CO flux, thus with the molecular gas 
content, it is noteworthy that the mean FIR surface density de¬ 
clines slightly from the Sa to the Sc type (Roberts & Haynes, 
1994). It should be therefore a fair approximation to assume 
that the relevant, HI 4- H 2 , gas content in the innermost, say 
0.5 kpc scale, region is not significantly (for our purpose) de¬ 
pendent on the morphological type. On the other hand, along 
the sequence Sa-Sc, the bulge component relative to the disk 
is largest in the earliest type and decreases thence. In addition, 
there exist a fairly close relationship between bulge luminosity 
and central BH mass (Ferrarese & Ford 2005). Thence, if the 
AGN were (statistically) correlated with Mbh, one should 
not expect the existence of the trend given by ^(Lx) 

This statement follows, of course, only if one refers to the 
gas content alone, and the gravitational influence of the BH 
(and the bulge) on the geometrical thickness of the gas within 


the first few hundred parsec from the center is disregarded. In 
the following, it is argued that this influence might indeed shed 
some light on the anticorrelation that one would like to explain. 
To illustrate the point, an “idealized” design is adopted, where 
the molecular gas is confined within a self-gravitating, starry 
galaxy disk, whose surface density E remains practically con¬ 
stant within a distance from the center of the order of one kpc. 
The scale height He of the molecular clouds is regulated by 
the disk gravity and their velocity dispersion. When closing in¬ 
wards, the bulge and the BH gravity adds to that of the starry 
disk, inducing a progressive shrinking of He- In geometrical 
terms this picture is quite similar to the one which applies to 
the physical models of the torus, summarized in the previous 
Section. In physical terms it is simpler, because outside a cer¬ 
tain radial distance the collisional and tidal effects have no se¬ 
rious impact and can be disregarded. 

The picture is characterized by the following parameters. 

- The number density no of the clouds at Z=0. 

- The scale heigth defined according to 

n(Z) = noe^*^' (3) 

- The clouds idealized as single spherical entities. Those 
which matter most for the absorption have the largest mass, 
i.e. the clouds that in our Galaxy dominate the mass func¬ 
tion (where most of the mass resides). Thus single val¬ 
ues for mass, radius, internal density are used: me, re, pe- 
Further parameters are the velocity dispersion (Je and the 
internal temperature Te- 

- The clouds are assumed to be in pressure equilibrium with 
an external gas (namely < mjeans), whose pe and Te are 

such that PeTe = PeTe, and He oc HeTe! (Je > He- 

Within at most 1 kpc both He and He are assumed constant 
and regulated by the disk force per unit mass, Kz- 

In this Section, the gravitational effects of a central BH (and 
the related Bulge mass) are described. Down to where, along R, 
it can be assumed (lacking obvious reasons to the contrary) that 
(Je and Te remain constant (as well as the Bulge density), this 
addition shall result into a change in He, namely this quantity 
will decrease with decreasing R, as soon as the vertical attrac¬ 
tion exerted by the BH will start dominating. The profile He(R) 
will go through an inflection at Rm/i, which regulates the cone 
opening angle, and increases with Mbh- Thus, if there were a 
statistical correlation between Mbh and Lx, the consequence 
would be a trend of the type given in Sect. 2, which is encour¬ 
aging. 

The force per unit mass exerted in the vertical, Z, direction 
by the matter in the disk can be expressed, within about Z=100 
pc, as 

Kz,d = -4nGpdZ (4) 

where pd is the approximately constant density of the matter in 
the disk, mainly provided by stars. 

The equivalent force exerted by the BH, in the approxima¬ 
tion Z <sc R, is given by 



( 5 ) 
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Similarly, the force exerted by the Bulge is 


Kz,b - -G 


Mb(R)^ 



(6) 


where Mh(R) is the bulge mass within the radius R. We have 
assumed a spherical symmetry, and a costant bulge density pi, 
within R (see the Appendix, where pt, is calculated as a function 
of the Black Hole mass). 

From the equation that links the clouds velocity dispersion 
and scale height to the sum of the three contributions to Kz 


of clouds along a line of sight. To this end, the vertical distri¬ 
bution of the clouds number density is assumed to follow an 
exponential law 

n{Z,R) = no{R)e^^T^^'' (10) 

where HdR) is given by eq. |3 and no is defined through the 
surface density 2, assumed constant with R: 

/^OO 

2 = 2mJ no{R)e~^T!^^^dZ (11) 

Jo 


pO 

cr^ - 2 I {Kz,d + Kz,bh + Kz,b)dZ 
Jhc 


where Me is the mass chosen to represent the single cloud. 
Thus: 


one infers 


Hc(R) = 


cr"- 


AnGpd + 


^TrGpb 




1/2 


( 8 ) 


3 


The distance of the point of inflection is then immediately 
derived to be 


Rinfl - 


1 


Mbh 


yAnpd+ 


1/3 


(9) 


niZ,R) = no-J^e ( 12 ) 

^C\^) 

where He is the scale height at » Rinfi- Except for the pa¬ 
rameter Z, all the other adopted quantities are given in Table 
1. In this table also a tidal radius is given, namely the distance 
within which the single cloud would be disrupted by the tidal 
force exerted by the BH, namely 


R tidal — 


/ 3Mbh 

ylnpcniH 


1/3 


(13) 


To be noted that at R=2Rinfi the gravitational pull of the 
BH in the Z direction is equal to that of the matter in the disk 
and the bulge. For values of R much lower than 2Rinfi, where 
the BH dominates, we shall disregard the contribution of the 
shrinking disk gravity (that is the increase of pd with decreas¬ 
ing R), and the profiles of He shown in Figure|2]correspond to 
a choice of pd equal to the value typical of a galaxy like our 
own, namely 0.85 solar masses per cubic parsec, a velocity dis¬ 
persion cr = 10 km i ', and four values of Mbh (and thence of 
Pb). The opening angle is very sensitive to the last parameter, 
and depends on the value of increasing from 12 pc when 
Mbh =10^ Mq to 226 pc when Mbh = IO^Mq. 


z 



Fig. 2. Scale height of the clouds as a function of R when 
Mbh=10® Mq, 10^ Mq, 10*^ Mo and 10® Mq. 


In order to evaluate the effective obscuration, this geometri¬ 
cal effect must be accompanied by assumptions on the number 


where pe is the internal density of the cloud, assumed equal 
to 10"^ hydrogen molecules cm“^. It should be noted that Rndai 
is of the same order of Rmfi for a BH mass of 10® Mq, and 
lower for higher masses. The adopted cloud size is such that the 
corresponding column density is equal to 10^^ cm“^. Because 
this is the value of N// chosen to discriminate between Seyferts 
with no absorption and Compton-thin absorbed ones, we have 
then simply to calculate the probability that one single cloud 
happens to lie along the line of sight. It is easy to demonstrate 
that it does not make a difference if such clouds, or smaller 
ones, where congregated in complexes rather than dispersed, 
as assumed, provided that the number of complexes remains 
large. 

The line of sight is defined through the angle 6 given by 

I = tgO, (14) 

and the number density as a function of 6 and R obtains from 

eq.El 

n(0,R) = no——e (15) 

For a given 9 the number of encountered clouds is given by 

pRo,,i 

N^nrl n(0,R)dR, (16) 

jRin 

where tt is the geometric cross section of a cloud, 
R„,i,=4Hc/tg6 and R,„ is consistently adopted equal to Rddai- 
Figure [^illustrates the number of clouds encountered as a 
function of 6, for four values of the Mbh and Z equal to 200 

Mq pc-2 
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disk matter density 

pd = O.S4MqPC * 

clouds velocity dispersion 

cr=10 km/s 

scale height of the clouds 

H=50 pc 

typical cloud inner density 

pc = 10'*cm“* 

typical cloud radius 

rc=0.167 pc 

typical cloud column density 

Nh,c = 

typical cloud mass 

Me = 4.35Mo 

typical cloud temperature 

r, = 10 K 

external gas temperature 

Tc = 10'' K 

external gas density 

pc = 0.1 cm“* 

cloud number density at z=0 

no = 'LIlHcMc 

tidal radius R for Mbh = 10^'’Mq 

R„rf„/=128 pc 

tidal radius R for Mbh = 10* Mq 

Rndai=60 pc 

tidal radius R for Mbh = 10'Mq 

Rtidai=2& pc 

tidal radius R for Mbh = 10'’Mq 

Rtidal=i3 pc 

inflection radius R for Mbh - IO^Mq 

Ri„f,=226 pc 

inflection radius R for Mbh = 10* Mq 

f^,>,//=101 pc 

inflection radius R for Mbh = 10'Mq 

Ri„ fi=39 pc 

inflection radius R for Mbh = 10^ Mq 

R/„//=12 pc 


Table 1. Adopted values for the parameters. See text for details. 



theta 


Fig. 3. Number of clouds along the line of sight as a function 
of 0 for four values of Mbh, when S=200 M© pc“^ 


By defining 0* as the angle where this number equals 1, it 
is straightforward to estimate the solid angle, hence the fraction 
of sources Ns where the line of sight is bound to cross a column 
of at least 10^^ cm“^, 


NsilogNn > 22) 

^s,tot 


cos(90 - 0*) 


(17) 


as a function of two parameters, the BH mass and the surface 
density 2 of the molecular gas within Rout- 

In Figure |4] (left panel), the fraction of absorbed AGN is 
shown for four different values of the BH mass (10^, 10^, 10^ 
and 10® Mo) and several choices of S. In order to draw a com¬ 
parison with the observational evidence, where the observable 


is the X-ray luminosity, it is necessary to adopt a relation be¬ 
tween Lx and Mbh- The luminosity-dependent bolometric cor¬ 
rection in Marconi et al. (2004; see their eq.(21) and Fig. 3) and 
a luminosity at 0.1 times the Eddington limit (see e.g. Peterson 
et al. 2004) have been adopted. The results are given in Figure|4] 
(right panel), where the best-fit relationship between the frac¬ 
tion of absorbed AGN and X-ray luminosity obtained by La 
Franca et al. (2005) is shown for comparison. 

The dependence of the slope of the anticorrelation on the 
adopted value of 2 is noteworthy. While for values of this pa¬ 
rameter larger than about 80-100 M© pc“^ the slope changes 
rather little, and remains gratifyingly close to the observed one, 
for smaller values it becomes shallower, and the more so the 
lower the value. This comes out simply because, in order to 
intercept a cloud, one needs to integrate well beyond the inflec¬ 
tion radius: in other words, despite the shrinking in He due to 
the black hole, the number density remains too small to have 
any sizeable effect within 2/?,>,//, and the required absorption is 
obtained only for line of sights grazing the plane of the disk. 

5. Discussion and conciusions 

The main result of this paper is that the anticorrelation be¬ 
tween ^(Lx) (namely the fraction of AGN with N^, as mea¬ 
sured by their X-ray spectra, greater than 10^^ cm“^) and Lx in 
Compton-thin objects is at least qualitatively reproduced pro¬ 
vided that; a) there is a sufficiently large number of molecular 
clouds within the radius where the BH gravitational influence 
is dominant; b) there is, in a statistical sense, a correlation be¬ 
tween the BH mass and the AGN luminosity. These two condi¬ 
tions are now briefly discussed. 

As shown in Fig. |3 the slope and the normalization of 
the anticorrelation are consistent with the observational results 
for values of S > 150-200 M© pc“^, if a bolometric lumi¬ 
nosity 0.1 times the Eddington luminosity is assumed (as it 
seems to be the case in the local universe, e.g. Peterson et al. 
2004). Thus, the first question arising is what are the typical 
values of this quantity in the innermost region of spiral galax¬ 
ies (the hosts of Seyfert active nuclei). The relevant radius of 
this region, 2Ri„fi, lies between about 25 and 450 pc (for Mbh 
from one and one thousand million solar masses), which corre¬ 
sponds to 0.5-9 arcsec at a distance of 10 Mpe, and requires the 
best angular resolution currently achievable, with interferome¬ 
try at the 3 mm CO 1=1-0 line, on relatively nearby galaxies. 
In this respect, the BIMA Survey of Nearby Galaxies (BIMA 
SONG) at Hat Creek, CA, with the 10-element Berkeley- 
Illinois-Maryland Association (BIMA) millimeter interferome¬ 
ter (Heifer et al. 2003) is the best available for our purpose, be¬ 
cause the sample was deliberatly selected without reference to 
CO or infrared brightness. It consists of 44 spiral galaxies of all 
morphological types (the Sa type are somewhat undersampled 
relative to the later types). The angular resolution of 6 arcsec 
corresponds to 360 pc (that is 180 pc radius for the central re¬ 
gion) at the average distance of the galaxies, 12 Mpc. Prom the 
CO surface brightness distribution (see Table 5 in Heifer et al. 
2003) it turns out that the peak brightness is attained within the 
central 6 arcsec in 20 out of 44 galaxies. In these 20 objects the 
corresponding surface density l,centr lies between 15 and 1854 
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Fig. 4. Expected behaviour of the fraction Nj(logNH > 22)/Ns^tot, for several values of E, as a function Mbh (left hand panel) and 
of L(2-10 keV) (right hand panel), with the latter derived assuming LBoi/LE(id=0.l (e.g. Peterson et al. 2004) and the luminosity- 
dependent bolometric correction in Marconi et al. (2004). In the right hand panel, crosses represent the observed fraction (cor¬ 
rected for selection effects, and adapted from La Franca et al. 2005) of obscured AGN in the local universe. 


Mo pc“^, with most of them (18) above 90, and the other two 
around 15 Mq pc“^. Notably, among the other 24 galaxies, 10 
have 1.centr > 50 Mo pc“^, with values up to 1151. The galaxies 
with Y.centr ^ 150 Mq pc“^ are 22/44, namely 50%. The straight 
mean value of this 22 galaxies is 511 Mq pc“^, notably close to 
the value of ~500 Mq pc“^ of the 300 pc molecular disk in our 
own galaxy (Glisten 1989). 

The sample contains 8 objects classified as Seyfert galaxies, 
6 of them with 1.centr greater than 90 and up to 466 Mq pc“^, 
1 (NGC 4725) with = 21 Mq pc'^, 1 (NGC3031) with 
'^centr below the detection limit; the proportion of galaxies with 
an active nucleus, which possess a high value of 'E.cemr, appears 
to be, within the small statistics, similar to that in the general 
population. 

The fraction of objects which meet our requirements is en¬ 
couraging, especially because it does not require a special cor¬ 
relation with the nuclear activity, and represents a relatively 
large scale (compared to the close circumnuclear environment) 
property fairly common in spiral galaxies. In other words, the 
condition required is already in place in about half of the galax¬ 
ies before a high rate nuclear activity lights up, which in turn 
is unlikely to significantly influence the preexisting condition. 
Much higher resolution and sensitivity surveys are needed to 
explore more distant objects, and therefore draw more solid 
conclusions on the validity of our hypothesis: to pursue this 
goal an array of millimeter telescopes like ALMA is necessary. 

Concerning the increase of the fraction of absorbed sources 
with redshift (La Franca et al. 2005), an interesting possibil¬ 
ity is that it might be due to an average increase of M with 
the redshift, if one were allowed to assume that the molecular 
content and distribution in spiral galaxies is independent of z. 
Indeed, McLure & Dunlop (2004) found that the Lboi /LEdd ra¬ 


tio raises from about 0.15 at z ~0.2 to 0.5 at z ~2. We intend 
to investigate this issue further, in particular within the frame 
of astrophysically self consistent models, which associate the 
growth of the supermassive BH in galaxies with the evolution 
in cosmic time of the AGN population (e.g. Menci et al. 2004). 
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where pi, is the central density and Rc is the core radius, if it is 
assumed that Rc is proportional to Ri,, it follows that Rc oc 
After normalization to the values which apply to the Milky 
Way, namely Rc,mw= 400 pc, M X lO^A/©, one obtains: 
Mbh = lO^Mo Rc = 180 pc 
Mbh = lO’Mo Rc = 700 pc 
Mbh = 10*^Mo Rc = 2.7 Kpc 
Mbh = 10®Mo 7?^ = H Kpc 

For these values of M.bh Rc turns out to be much larger than 
two times dehned in Sect. 4 and given in Table 1. Thus 
it can be concluded that the stellar bulge density within 
remains practically constant. Its value, in the above approxi¬ 
mations, is proportional to M^J^^, and we calculated it for each 
BH mass by simply rescaling the Milky Way value: 

Mbh = lO^Mo pt = 15.1 Mo/pc^ 

Mbh = lO^Mo pb = 2.4 MJpc^ 

Mbh = 10*^Mo pb = 0.4 Mo/pc^ 

Mbh = IO^Mq pb = 0.06 M^lpc^ 


Appendix A: Effects of the stellar bulge 

The BH mass is known to correlate with the velocity disper¬ 
sion, cTb, in the bulge: 


Mbh 

lO^Mo 


(1.66 + 0.24) 


( r 


with a - 4.86 ± 0.43 (Ferrarese & Ford 2005). 
With respect to the virial mass of the bulge 

GMb (T^Rb, 


(A.l) 


(A.2) 


the BH mass turns out (Merritt & Ferrarese 2001) proportional 
to Mft, Mbh - lO^^Mfe. From ea. lA.ll and lA.2l one immediately 
obtains 


Mb Mbh .A-i/a 

^BH 


^BH 


(A.3) 


or, rounding off a to 5, 

R bulge ^ K^BH' 


(A.4) 


Adopting for the density prohle of the bulge the relation¬ 
ship: 


p(R)^ 


Pb 

1 + (RIRc? 


(A.5) 






